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ZEUS Results on the Measurement and 
Phenomenology of F<i at Low x and Low Q 2 



Abstract 

Measurements of the proton structure function F 2 for 0.6 < Q 2 < 17 GeV 2 and 
1.2 x 10~ 5 < x < 1.9 x 1CT 3 from ZEUS 1995 shifted vertex data are presented. 
From ZEUS F 2 data the slopes dF 2 /d\nQ 2 at fixed x and din F 2 / dln(l / x) for 
a; < 0.01 at fixed Q 2 are derived. For the latter E665 data are also used. The 
transition region in Q 2 is explored using the simplest non-perturbative models and 
NLO QCD. The data at very low Q 2 < 0.65 GeV 2 are described successfully by 
a combination of generalised vector meson dominance and Regge theory. From a 
NLO QCD fit to ZEUS data the gluon density in the proton is extracted in the 
range 3 x 10~ 5 < x < 0.7. Data from NMC and BCDMS constrain the fit at large 
x. Assuming the NLO QCD description to be valid down to Q 2 ~ 1 GeV 2 , it is 
found that the qq sea distribution is still rising at small x and the lowest Q 2 values 
whereas the gluon distribution is strongly suppressed. 
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1 Introduction 



The rapid rise of the proton structure function F2 at low x, measured at HERA M, 
continues to generate a lot of interest. In particular the persistence of the strong rise to 
small values of Q 2 and the apparent success of the perturbative QCD (pQCD) descrip- 
tion of the data down to Q 2 values approaching 1 GeV 2 raise new challenges for our 
understanding of QCD. HERA also allows study of the 'transition region' as Q 2 — » in 
which pQCD must break down. The theoretical context for our study of pQCD and the 
transition region is outlined in Sec. ^. 

With data taken during the 1995 HERA run the ZEUS experiment has achieved a 
significant increase in the kinematic coverage for low x and low Q 2 inelastic neutral 
current positron-proton scattering. The coverage for Q 2 between 0.11 and 0.65 GeV 2 
was made possible with the installation of a small electromagnetic sampling calorimeter, 
the Beam Pipe Calorimeter (BPC), at small positron scattering angles and results on 
the proton structure function F2 and the total j*p cross-section have been published [|j . 
In Sec. H of this paper we report on further measurements of F% in the Q 2 range between 
0.6 and 17 GeV 2 . The data were obtained from runs in which the interaction point was 
shifted away from the main rear calorimeter thus extending its small-angle coverage for 
scattered positrons. These data fill the gap in Q 2 between the BPC and the 1994 ZEUS 
measurements ||, Q. Taking all three data sets together, the ZEUS experiment has 
measured F2 over the kinematic region 0.1 < Q 2 < 5000 GeV 2 , 2 x 10~ 6 < x < 0.5. The 
coverage of the kinematic plane by the ZEUS data sets is shown in Fig. |l[ 
The very low Q 2 data are discussed in Sec. || using generalised vector dominance and 
Regge theory and it is established that the ZEUS data with Q 2 < 0.65 GeV 2 are 
well described by such approaches. In Sec. || the slopes dF^jd In Q 2 at fixed x and 
d\a F2 / d\n.(l / x) at fixed Q 2 are derived from the combined ZEUS data sets. In Sec. ^ 
the ZEUS F2 data, together with fixed target data at large x, are fit using next to lead- 
ing order (NLO) QCD to determine the gluon momentum density. The increased range 
and precision of the ZEUS F2 data allow a more precise extraction of the gluon density 
at low x compared to our earlier results using the 1993 ZEUS data §. In Sec. the 
properties of the pQCD description of the F2 and slopes data at low Q 2 are explored 
in more detail. Our conclusions are summarised in Sec. |8[ Tables containing F2 values 
and other data are given in the Appendix. 



2 Phenomenology of the low Q 2 region 

We use NLO pQCD and the simplest non-perturbative models to explore the transition 
region in Q 2 . The standard NLO DGLAP equations || give the Q 2 evolution of parton 
densities, but do not prescribe their functional form in x Q at the starting scale Qq. At 
Qq, the small x behaviour of parton momentum densities xf(x) may be characterised by 
the exponent 5 where xf(x) ~ Ax 6 . For 5 > a parton momentum density either tends 
to zero or is constant as x — ► (non-singular) while for 5 < a parton density increases 
as x — > (singular). One way to understand the rise of F2 at low x is advocated by 
Gliick, Reya and Vogt (GRV) |7|, || who argue that the starting scale for the evolution 
of the parton densities should be very low (0.34 GeV 2 ) and at the starting scale the 



1 Bjorken x and the negative squared momentum transfer, Q 2 , are denned in Sec. 3.2 
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parton density functions should be non-singular. For Q 2 > 1 GeV 2 , the observed rise 
in F2 is then generated dynamically through the DGLAP evolution equations. On the 
other hand, at low x one might expect that the DGLAP equations break down because 
of large ln(l/x) terms that are not included. Such terms are taken into account by the 
BFKL formalism |J, which in leading order predicts a rising F2 at low x. The rise comes 
from a singular gluon density with 5 in the range —0.3 to —0.5. Recent work on BFKL 



at NLO has shown that the corrections to the LO value for 5 are large [1C] and reduce 
the predicted rise in F%, though quite how large the reduction should be is still under 
discussion Clearly accurate experimental results on F2 and 5 at low x are of great 
interest. More details on the many alternative pQCD approaches to the low x region 
may be found in [fi^, O. 



As Q 2 decreases as increases and pQCD will eventually break down. Then non- 
perturbative models must be used to describe the data. At low x the lifetime of the 
virtual photon in the proton rest frame is large compared to the "f*p interaction time 



[14]. Inelastic e p scattering may then be viewed as j*p scattering, with the total 7*p 



cross-section given by (] 

<?L P (W 2 , Q 2 ) =*t + *l* ^F 2 {x, Q 2 ), (1) 

where W ~ y/Q' 2 /x is the centre-of-mass energy of the j*p system and ot and gl 
are the cross-sections for transversely and longitudinally polarised virtual photons re- 
spectively. We consider two non-perturbative approaches, the vector meson dominance 
model (VMD) and Regge theory. 

VMD relates the hadronic interactions of the photon to a sum over interactions of the 
P°,uj and <p vector meson states |l5|, [Tq ]. To accommodate deep inelastic scattering data 
the sum has to be extended to an infinite number of vector mesons giving the generalised 



vector dominance model (GVMD) [17]. Following the assumptions in [17], axiW 2 ,Q 2 ^ 



is related to a' yp (W ), the total photoproduction cross-section by 

r c M 2 , ^ r v M* 



a T {W z ,Q 2 ) 



(M 2 + Q 2 ) + , (M£ + Q 2 ) 2 



a™(W 2 ), (2) 



where Mq is the lower cutoff of the continuum vector states and rc,ry are constants 
satisfying the normalisation condition tq + XV r V = 1 a t Q 2 = 0. A similar expression 
to Eq. @ may be written for <jl{W 2 , Q 2 ), but with additional Q 2 dependence to ensure 
that it vanishes as Q 2 — > ]17|]. The GVMD approach has recently been revived in the 
context of low Q 2 HERA data by Schildknecht and Spiesberger |l8[| . We use a simplified 
form of GVMD to study the consistency of the data in the ZEUS BPC region and its 
extrapolation to Q 2 = 0. 



Regge theory [pLS| , 20] provides a framework in which the energy dependence of hadronic 
total cross-sections is of the form a ~ ^ r (3 r s ar ~ l where yfs is the centre-of-mass energy, 
ol t the intercept of the Regge trajectory and (3 r a process dependent constant. The a r 
are universal and can in principle be determined from the spectrum of meson states. 
However, for the dominant trajectory describing total cross-sections at high energies 



2 Considering virtual photon exchange only. Since we are working at small x, terms depending on 
MfjX 2 have been ignored. 
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(known as the Pomeron), this has not yet been possible. The Pomeron intercept, ap, 
is determined by fitting high energy total cross-section data. Donnachie and Landshoff 
(DL) |2l|] have used a two component Pomeron+Reggeon approach to give a good de- 
scription of hadron-hadron and photoproduction total cross-section data over a wide 
range of energies with ap of about 1.08. They have extended their approach to 7*p 
total cross-sections [|^] by keeping the Regge intercept independent of Q 2 but assuming 
a simple Q 2 dependence for the coupling term, which becomes P r m 2 /(m 2 + Q 2 ) where 
m 2 is again determined by fitting to data. 

Neither the non-perturbative VMD and DL approaches nor pQCD can be expected to 
describe the Q 2 behaviour of F2 over the complete range from photoproduction to very 
large Q 2 deep inelastic scattering (DIS). Many models combining various aspects of 



these approaches have been applied to the data in the transition region [13, 23, A 



comparison of the ZEUS BPC data with some of the models was given in We use 
the DL Regge model ^ and the 1994 parton densities (GRV94) of the GRV group 
H as 'benchmarks' in this paper as they were available before the recent precision 
measurements were made. 



3 Measurement of F2 with shifted vertex data 

The shifted vertex data correspond to an integrated luminosity of 236 nb _1 taken in a 
special running period in which the nominal interaction point was offset in the proton 
beam direction by +70 cm, i.e. away from the main rear calorimeter (RCAL). The 
measurement follows previous analyses described in more detail in |j. However, 
compared to the earlier shifted vertex analysis [||, for the 1995 data taking period 
the RCAL modules above and below the beam were moved closer to the beam, thus 
extending the shifted vertex Q 2 range down to 0.6 GeV 2 . The basic detector components 
used are the compensating uranium calorimeter (CAL), which has an energy resolution, 
as measured in the test beam, of a/E = 18%/-\/ E(GeV) for electromagnetic particles 
and a/E = 35%/-*/ -©(GeV) for hadronic particles. The tracking chamber system is used 
to determine the position of the event vertex. The small angle rear tracking detector 
(SRTD) consists of horizontal and vertical scintillator strips covering the region around 
the RCAL beam hole, i.e. the region of positron scattering angles of low Q 2 events. It is 
also used as a preshower detector for the RCAL and has a position resolution of 0.3 cm. 
The luminosity is determined from the positron-proton bremsstrahlung ep — > ep7 where 
the radiated photon is measured in a lead-scintillator calorimeter (LUMI) positioned 
at Z = —107 m. There is an associated electron calorimeter (LUMI-E), positioned at 
Z = —35 m, which is used for tagging photoproduction events. The uncertainty in the 
luminosity measurement is 1%. 



3.1 Monte Carlo simulation 

Monte Carlo (MC) events are used to correct for the detector acceptance, resolution and 



the effect of initial state radiation. In the framework of DJANGO [25] the generator 



HERACLES [26] is used to simulate neutral current DIS events including first order 

3 The ZEUS coordinate system is denned as right handed with the Z axis pointing in the proton 
beam direction, and the X axis horizontal, pointing towards the centre of HERA. The origin is at the 
unshifted interaction point. 
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electroweak radiative effects. The hadronic final state is simulated using the ARIADNE 
p7| program which implements the colour-dipole model. A parameterisation of the F2 
structure functions based on the results published in § is used and Fl is set to zero. 
The MC event sample is 1.6 times that of the data and the events are passed through the 
same offline reconstruction software as the data. Simulated photoproduction background 
events are generated using the program PYTHIA [28] with a cross-section given by the 
ALLM [29] parameterisation. 



3.2 Kinematic reconstruction 

The reaction e + (k) + p(P) — ► e + (k') + X at fixed squared centre-of-mass energy, s = 
(k + P) 2 , is described in terms of Q 2 = —q 2 = —{k — k') 2 and Bjorken x = Q 2 /(2P ■ q). 
At HERA s « 4E e E p , where E e = 27.5 GeV and E p = 820 GeV denote the positron 
and proton beam energies. The fractional energy transferred to the proton in its rest 
frame is y = Q 2 /(sx). 

The kinematic variables are reconstructed from the measured energy, E' e , and scattering 
angle, 9 e , of the positron (the 'electron method'), 

E' 

Q 2 = 2E e E' e (1 + cos e ) and y = l §• (1 - cos B e ) . 

2E e 

This method gives the best resolution in the region of interest at high y and low Q 2 . 
Scattered positrons are identified by a neural network based algorithm [Q], with an 
efficiency of about 90% at positron energies of 10 GeV, increasing to 100% at 20 GeV. 
The measured energy is corrected for energy loss in inactive material in front of the CAL 
using the signals in the SRTD scintillators. The uncertainty in the measured energy is 
estimated to be 2% at 10 GeV decreasing to 1% at 27.5 GeV. The positron impact 
position on the CAL measured with the SRTD together with the event vertex position 
measured with tracks in the CTD gives the positron scattering angle 8 e . For events 
outside the fiducial volume of the SRTD, the CAL position determination is used. For 
events in which the event vertex cannot be reconstructed, the vertex is set to the mean 
value of the vertex distribution. 

The variable yjB = J2{Ei(l — cos 6i)/(2E e ), where the sum runs over all CAL cells 
except those belonging to the scattered positron, gives a measurement of y with a good 
resolution at low y. A cut ?/jb > 0.04 is imposed to limit event migrations from low y, 
where the resolution of the electron method is poor, into the bins at higher y. 



3.3 Event selection 

Data are selected online by a three level trigger system. At the first level a certain energy 
deposit in the CAL is required and cuts on the arrival times of particles measured in 
the SRTD are imposed. At the second level the condition 5 = J2i ^t(l ~~ cos^) > 
29 GeV — 2E^ has to be fulfilled, where the sum goes over all CAL cells with energies 
Ei and polar angles 9i, and E^ is the energy measured in the LUMI detector. This 
cut significantly reduces the photoproduction background as 5 = 2E e (55 GeV) for a 
fully contained DIS event. At the third level a full event reconstruction is performed. 
A reconstructed positron with an energy greater than 4 GeV and a CAL impact point 
outside a box of 24(A) cm x 12(Y) cm centered on the RCAL beam hole is required. 
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Also, S has to be greater than 30 GeV and event times measured in the CAL are required 
to be consistent with an e + p interaction at the nominal shifted interaction point. 
The offline event selection cuts are: 

• Positron finding as described above, including the requirements on the impact 
point and on ?/jb > 0.04, with a corrected positron energy E' e > 10 GeV. This 
ensures a high efficiency for positron finding and removes events at very high y, 
which suffer from large photoproduction backgrounds. 

• The positron impact point on the CAL is required to be outside a box of 26 cm x 
16 cm around the RCAL beam pipe hole to ensure full shower containment in the 
CAL. 

• 35 GeV < 5 < 65 GeV, in order to further reduce photoproduction and beam-gas 
related backgrounds. This cut also removes events with hard initial state radiation. 

• For events with a tracking vertex, the reconstructed Z coordinate of the vertex is 
required to lie within 40 cm < Z vertcx < 160 cm. The acceptance is extended to 
large Z values to accommodate events from satellite bunches, i.e. proton bunches 
that are shifted by 4.8 ns with respect to the primary bunch crossing time, resulting 
in a fraction of ep interactions occurring displaced by an additional 72 cm. 

A total of 62000 events pass the cuts. 



3.4 Background estimation 

The background from beam-gas interactions is about 1% as determined from unpaired 
positron and proton bunches. 

The main background comes from photoproduction events, where the positron escapes 
along the beam line and a mis-identified positron (mainly electromagnetic showers from 
7T° decays) is reconstructed in the CAL. The amount of this background is determined 
using the MC simulated photoproduction event sample. In total it is a small effect which 
is only significant at small values of E' e and Q 2 as shown in plots (d) and (e) of Fig. ||. 
In a small fraction of real photoproduction events the positron is detected inside the 
limited acceptance of the LUMI-E electron tagger. These events are used to cross check 
the MC background estimate. Both results agree within 20%. 



3.5 Determination of F 2 



In the Q 2 range of this analysis the double differential cross-section for single virtual- 
photon exchange in DIS is given by 



d 2 a 



2vra 2 



dxdQ 2 xQ 4 



2(l-l/) + 



V 



l + R 



F 2 (x,Q 2 ) l + 5 r (x,Q 2 



(3) 



where R is related to the longitudinal structure function Fl by R = Fi,l{F 2 — Fl) and 
5 r gives the radiative correction to the Born cross-section. For the kinematic range of 



this analysis 8 r is at most 10%. For R we take values given by the BKS model [30]. An 
iterative procedure is used to extract the structure function F 2 . Data and MC events 
are binned in the variables y and Q 2 . In a bin- by-bin unfolding procedure the MC 
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differential cross-section is adjusted to describe the data using a smooth function for F<i- 
The re-weighted MC events are then used to unfold Fi again, until after 3 iterations 
the changes to F% are below 0.5%. The statistical errors of the F2 values are calculated 
from the number of events measured in a bin and the statistical error on the acceptance 
calculation from the MC simulation. 

The quality of the description of the data by the re-weighted MC is shown in Fig. ||, 
which displays distributions of the following quantities: (a) 5 or E — Pz as defined in 
Sec. (b) the Z-position of the primary vertex; (c) the positron scattering angle 9 e ; 
(d) the energy, E' e , of the scattered positron; (e) log 10 Q 2 ; (f) log 10 y. The agreement 
between data (filled circles) and simulated data (open histograms, normalised to the 
luminosity of the data) is generally good. 

3.6 Systematic uncertainties 

The systematic uncertainties of the measured F% values are determined by changing 
the selection cuts or the analysis procedure in turn and repeating the extraction of F^. 
Positive and negative differences, AF2, are added in quadrature separately to obtain the 
total positive and negative systematic errors. The systematic checks and errors are: 

• A shift of the horizontal and vertical position of the SRTD by ±0.5 mm results in 
AF2/F2 of at most » 2%. 

• The two halves of the SRTD are moved with respect to each other in the horizontal 
and vertical directions by 1.0 mm, which results in AF2/F2 of ~ 4%. 

• The uncertainty in the positron energy calibration is estimated to be 2% at 10 GeV 
decreasing to 1% at 27.5 GeV giving a maximum AF2/F2 of 8%. 

• The hadron energy scale is uncertain to ±3%, causing a maximum AF2/F2 of 2%. 

• The uncertainty in the positron finder efficiency is estimated to be 2.5% at 10 GeV 
decreasing to 1% at 27.5 GeV, which gives AF2/F2 of at most 2%. 

• The number of MC events in satellite bunches is increased by 100% and decreased 
to 50%. This leads to a AF 2 /F 2 less than 2%. 

• The uncertainty in the determination of the photoproduction background by MC is 
estimated to be 30%. The resulting maximum AF2/F2 is about 7% in the highest 
y bins. 

• The vertex finding efficiency is between 75% and 95% depending on the kinematic 
region. To study the effect of differences in the efficiency between MC and data 
all vertices are assigned to the nominal interaction point. AF2/F2 is at most 6%. 

• A variation in the box cut from 26 cm x 16 cm to 25.6 cm x 15.6 cm or 26.4 cm x 
8.4 cm, giving a maximum change AF2/F2 of 7%. 

The acceptance for DIS events with a rapidity gap at low Q 2 is somewhat different from 
that of non-diffractive events due to the different energy flow. To check the effect on F2, 
the acceptance for diffractive events is first calculated using a separate diffractive MC.Q 

4 A modification of the ARIADNE MC adjusted to generate rapidity gap events as described in [0. 
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Using this and the measured fraction of rapidity gap events in each bin, the acceptance 
function is recalculated. The largest change to F2 is 2.5%. The data are corrected for 
this effect. Half the correction value is taken as the estimate of the systematic error, 
reflecting mainly the uncertainty in the fraction of diffractive events. 
In addition there is an overall normalisation uncertainty of 1.5%, due to the 1% error 
in the luminosity measurement and a 1% uncertainty in the trigger efficiency, which is 
not included in the point to point systematic error estimate. 

The F2 data cover the x range 1.2 x 1CT 5 - 1.9 x 1CT 3 in 12 bins of Q 2 between 0.6 
and 17 GeV 2 (ZEUS SVX95). The values for F2 and their systematic errors are given in 
Table |l] of the Appendix. Fig. [3| shows the results for F2 as a function of x in the 12 Q 2 
bins. In the lowest Q 2 bin data from ZEUS F% measurements at very low Q 2 using the 
BPC (BPC95) |2j are shown and at larger Q 2 those from the ZEUS94 measurements 
Also shown are data from the shifted vertex measurements by HI (HI SVX95) |3l| and 



fixed target data from E665 32]. There is good agreement between the different ZEUS 
data sets and between ZEUS, E665 and HI data in the regions of overlap. We note that 
the steep increase of F2 at low x observed in the higher Q 2 bins softens at the lowest 
Q 2 values. The curves shown will be discussed later. 



4 The low Q region 

We first give an overview of the Q 2 < 4.55 GeV 2 region. Fig. I shows the ZEUS 
cross-section data versus W 2 derived from the SVX95, BPC95 and ZEUS94 F 2 values 
using Eq. (|l|). Also shown are data from HI SVX95 and measurements of the total 
cross-section for scattering of real photons on protons at fixed target [33] and HERA 



energies [34]. The two curves shown are the predictions of the DL Regge model [22] 
and cr^ ot p calculated from the NLO QCD parton distributions of GRV94 ||. The DL 
model predicts that the cross-section rises slowly with energy oc W 2X , A = 1 — ap ~ 0.08 
and this behaviour seems to be followed by the data at very low Q 2 values, although 
the normalisation of the DL model is low compared to the ZEUS BPC95 data. Above 
Q 2 = 0.65 GeV 2 , the DL model predicts a shallower rise of the cross-section than the 
data exhibit. For Q 2 values around 1 GeV 2 and above, the GRV94 curves describe 
the qualitative behaviour of the data, namely the increasing rise of o"7 ot P with W 2 , as 
Q 2 increases. This suggests that the pQCD calculations can account for a significant 
fraction of the cross-section at the larger Q 2 values. 

For the remainder of this section we concentrate on non-perturbative descriptions of the 
ZEUS BPC95 data (0.11 < Q 2 < 0.65 GeV 2 ). Since the BPC data are binned in Q 2 and 
y we first rewrite the double differential cross-section of Eq. (|3|) (dropping the radiative 
correction factor) as 

Q 2 

where ol = Fl and gt + &l has been defined by Eq. (pi). The virtual photons 
4ir z a 

have flux T = a(l + (1 - y) 2 )/(27rQ 2 y) and polarisation e = 2(1 -y)/{l + (1 -y) 2 ). For 
the BPC data e lies in the range 0.55 — 0.99 but as s is fixed at HERA, e cannot be 
varied independently of x and Q 2 . Thus the experimently determined quantity is the 
combination gt + cgl- For simplicity we keep only the continuum term in the GVMD 
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expression of Eq. (g). At a fixed W the longitudinal and transverse 7*1? cross-sections 
are then related to the corresponding cross-section a2 p at Q 2 = by 



vl(W 2 ,Q 2 ) = £ 



M 2 M 2 + Q 2 Mq 



Q 2 Mq M 2 + Q 2 



M 2 



MW 2 ,Q 2 ) = j^fQ- 2 <(w% (5) 

where the parameter £ is the ratio a^J ' /oj, p for vector meson (V) proton scattering 
and Mq is the effective vector meson mass. Neither £ nor Mq are given by the model 
and they are usually determined from data. We set cjl to zero because £ is expected 
to be less than one [] and the factor in the square bracket in the expression for er^ is 
small (for Q 2 < 0.65 GeV 2 and Mq m p , it is less than 0.2). The Q 2 dependence 
of the BPC data, in 8 bins of W between 104 and 251 GeV, is fit with a single mass 
parameter Mq. The cross-sections, aQ P (W 2 ), are also fit at each W giving a total of 
9 parameters. The fit is reasonable (x 2 /ndf = 38.7/(34 — 9) = 1.55, statistical errors 
only) as shown in the upper plot of Fig. |5|. To estimate the systematic errors, the fit 
is first repeated for each systematic check on the BPC data, using data and statistical 
errors corresponding to that change. The systematic errors are then determined by 
adding in quadrature the changes from the nominal values of the parameters. As a final 
check on the stability of the results, the fit is repeated including the longitudinal term 
with £ = 0.4. The resulting changes in the values for the cross-sections are less than 
their statistical errors (more details are given in |3q1 ). The value obtained for Mq is 
0.53 ± 0.04(stat) ± 0.09(sys) GeV 2 . The resulting extrapolated values of cjq P are given 
in Table || of the Appendix and shown as a function of W 2 in the lower plot of Fig. [|, 
along with measurements from HERA and lower energy experiments. The BPC &q P 
values lie somewhat above the direct measurements from HERA. They are also above 
the prediction of Donnachie and Landshoff. It should be clearly understood that the 
<7q P values derived from the BPC are not a measurement of the total photoproduction 
cross-section but the result of a phenomenologically motivated extrapolation. 
The simple GVMD approach just described gives a concise account of the Q 2 dependence 
of the BPC data. To describe the energy dependence of the data we use a two component 
Regge model 

aZ(W 2 ) = AniW 2 )**- 1 + A P {W 2 ) a r- 1 

where P and R denote the Pomeron and Reggeon contributions. The Reggeon in- 
tercept olr is fixed to the value 0.5 which is compatible with the original DL value 



TJ, and recent estimates [p7|, 38 1. Fitting only the Pomeron term to the extrapo- 



lated BPC data (tr^) gives a P = 1.141 ± 0.020(stat)±0.044(sys). Fitting both terms 
to the real photoproduction data (with W 2 > 3 GeV 2 ) and BPC <7q P data yields 
ap = 1.101 =t 0.002(stat)±0.012(sys). Including, in addition, the two direct measure- 
ments from HERA @ gives a P = 1.100 ± 0.002(stat)±0.012(sys). At HERA energies 
the contribution of the Reggeon term is negligible. The values of ap are compatible 
with the DL value of 1.08 and the recent best estimate of 1.0964^Q' Qg 4 by Cudell et al. 



Prom studies of diffractive vector meson production and data on R(= ol/ot) in inelastic ep scat- 
tering £ is in the range 0.2 — 0.4 [Q. More recently ZEUS |3^] has measured ep epp at low Q 2 and 
found that R is about 0.4 at Q 2 = 0.5 GeV 2 . 
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The final step in the analysis of the BPC data is to combine the Q 2 dependence from 
the GVMD fit with the energy dependence from the Regge model 

4o?{W\Q 2 )= (j^Q^j (ArQV*)**- 1 + A P (W 2 r^). (6) 

The parameter Mq is fixed to its value of 0.53 found above and an is also kept fixed at 
0.5 as before. The 3 remaining parameters are determined by fitting to photoproduction 
data (with W 2 > 3 GeV 2 , but without the two original HERA measurements ) and the 
measured BPC data. We find A R = 145.0±2.0//b, A P = 63.5±0.9/xb and a P = 1.097± 
0.002 with x 2 /ndf = 1.12 (statistical errors only). If the two HERA photoproduction 
measurements are included the parameters do not change within their errors. The 
description of the low Q 2 F 2 data given by this model (ZEUSREGGE) is shown in 
Fig. H Data in the BPC region Q 2 < 0.65 GeV 2 are well described. At larger Q 2 values 
the curves fall below the data. Including ZEUS SVX95 data at successively larger Q 2 
values, we find that by Q 2 = 3.5 GeV 2 the x 2 /ndf has increased to 1.7. Also shown in 
Fig. g for Q 2 > 0.9 GeV 2 , are the results of a NLO QCD fit (ZEUSQCD) described in 
Sec. § 

To summarise, we have shown that the Q 2 dependence of the ZEUS BPC95 data at very 
low Q 2 can be described by a simple GVMD form. The resulting values of <Tq P , the cross- 
sections extrapolated to Q 2 = 0, are somewhat larger than the direct measurements at 
HERA. A two component Regge model gives a good description of the W 2 dependence 
of the data, with a Pomeron intercept compatible with that determined from hadron- 
hadron data. 



5 F2 slopes 

To quantify the low x behaviour of F 2 as a function of x and Q 2 , we calculate the slopes 
dlnF 2 /dln(l/x) at fixed Q 2 and dF 2 /dlnQ 2 at fixed x from the ZEUS SVX95, BPC95 



and ZEUS94 data sets. We use the ALLM97 parameterisation 3£] for bin-centering F 2 
data when necessary. 

5.1 The slope dlnF 2 /rfln(l/x) 

It is seen from Fig. [6| that the x slope of F 2 is small for small Q 2 and then starts to 
increase as Q 2 increases. At a fixed value of Q 2 and at small x the behaviour of F 2 
can be characterised by F 2 oc x~* e ff (giving \ e ff = d In F 2 /dln(l/x)), with \ e ff taking 
rather different values in the Regge and LO BFKL approaches. The value of X e ft as 



an observable at small x has been discussed by Navelet et al. |4(], 41 1 and data on X e ff 



with x < 0.1 have been presented by HI [42]. 



Using statistical errors only, we fit F 2 data at fixed Q 2 and x < 0.01 to the form Ax~ Xe ff . 
Refering to Fig. ffl, we are measuring X e jf from horizontal slices of F 2 data between the 
y = 1 HERA kinematic limit and the fixed cut of x < 0.01. As the x range of the ZEUS 
BPC95 data is restricted we include data from E665 |32[| . In each Q 2 bin the average 
value of x, (x), is calculated from the mean value of ln(l/x) weighted by the statistical 
errors of the corresponding F 2 values in that bin. For the estimation of the systematic 
errors, it is assumed that the systematic error analyses for each of the four data sets used 
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are independent. For a particular data set and a given systematic check, the F2 points 
in each Q 2 bin are moved up and down by the respective systematic error estimates and 
the fits repeated, keeping all other data sets fixed at their nominal values. The positive 
and negative shifts with respect to the central values of X e jf are added separately in 
quadrature to give the positive and negative systematic errors. 

Fig. |j] shows the measured values of \ e ff as a function of Q 2 , and the data are given 
in Table ^ of the Appendix. At very low Q 2 the errors on \ e ff are large because 
this region is below the lower limit of E665 data (see Fig. |). At Q 2 > 100 GeV 2 
the statistical error dominates. From the Regge approach of the previous section one 
would expect A e // ~ 0.1, independently of Q 2 . Data for Q 2 < 1 GeV 2 are consistent 
with this expectation. The points labelled DL, linked by a dashed line, are calculated 



from the Donnachie-Landshoff prediction 22] and as expected from the discussion of 



the previous section are somewhat below the data. The variation of the DL points with 
Q 2 is a consequence of averaging the model in a Q 2 bin over a variable range of x and 
hence W 2 . For Q 2 > 1 GeV 2 , A e // increases to around 0.3 at Q 2 values of 40 GeV 2 . 
Qualitatively the tendency of X e ff to increase with Q 2 is described by a number of pQCD 
approaches The points labelled GRV94, linked by a dashed line, are calculated from 
the NLO QCD GRV94 fit. Although the GRV94 prediction follows the trend of the data 
it tends to lie above the data, particularly in the Q range 3-20 GeV 2 . We shall return 
to this point later. For the predictions shown in Fig. [7] the same F% error weighted 
average in x at a given Q 2 is used as for the data. 

5.2 The slope dF 2 /d\nQ 2 

In QCD the scaling violations of F2 are caused by gluon bremsstrahlung from quarks 
and quark pair creation from gluons. In the low x domain accessible at HERA the 
latter process dominates the scaling violations. F2 is then largely determined by the 
sea quarks F2 ~ xS, whereas the di^/cHnQ 2 is dominated by the convolution of the 
splitting function P qg and the gluon density: cii^/d In Q 2 ~ ctsP qg <S> xg. This has been 
used by Prytz to relate xg directly to the measured values of di^/tilnQ 2 0. The 
importance of dF2/d\nQ 2 as a tool for studying the low x region was pointed out by 
B artels et al. Q. 



In order to study the scaling violations of F2 in more detail the logarithmic slope 
dF2/ 'dlnQ 2 is derived from the data by fitting F% = a + blnQ 2 in bins of fixed x, using 
only statistical errors. The ZEUS data sets used are the BPC95, SVX95 and ZEUS94. 
For compatibility with our NLO QCD fit a cut of W 2 > 10 GeV 2 is applied to the data. 
The F2 data are shown in bins of x as functions of Q 2 in Fig. [8|. The fits F2 = a + b In Q 2 
are also shown and the change of slope as x changes is visible from the plots. In each x 
bin the average value of Q 2 , (Q 2 ), is derived from the F2 statistical error weighted mean 
value of In Q 2 in that bin. Systematic errors on dF2 /d\nQ 2 are estimated following the 
procedure outlined in the previous section. The results for dF2/dlnQ 2 as a function of 
x are shown in Fig. |9] and are given in Table |4] of the Appendix. The differences in the 
sizes of the errors on di^/cHnQ 2 partially reflect the variations in Q 2 range as x varies 
(see Fig. |l| - we are taking vertical slices of the data to determine di^/cHnQ 2 )- For 
values of x down to 3 x 10 -4 , the slopes are increasing as x decreases. At lower values 
of x and Q 2 , the slopes decreases. If di^/dlnQ 2 values are plotted for fixed target data 
at similar values of Q 2 , the 'turn over' is not seen, but the data are at larger values of 
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x [45, 46]. The points linked by the dashed lines are again from the DL Regge model 
and the GRV94 QCD fit, in both cases calculated using the same F^ error weighted Q 2 
averaging as for the data. The failure of DL is in line with our earlier discussion but 
the fact that GRV94 does not follow the trend of the data when it turns over is perhaps 
more surprising. Naively it appears that the GRV94 gluon density is too large even at 
Q 2 values around 5 GeV 2 . We shall return to this discussion after we have presented 
the ZEUS NLO QCD fit to which we now turn. 



6 NLO QCD fit to F2 data and extraction of the gluon 
momentum density 

In this section we present a NLO QCD fit to the ZEUS94 data Q and the SVX95 data of 
this paper. We do not attempt a global fit to parton densities, but concentrate on what 
ZEUS data allow us to conclude about the behaviour of the gluon and quark densities 
at low x. 

To constrain the fits at high x, proton and deuteron F% structure function data from 
NMC @ and BCDMS @ are included.^] The following cuts are made on the ZEUS and 



the fixed target data: (i) W 2 > 10 GeV 2 to reduce the sensitivity to target mass |49 l 
and higher twist j5(| contributions which become important at high x and low Q 2 ; 
(ii) discard deuteron data with x > 0.7 to eliminate possible contributions from Fermi 



motion in deuterium [51]. The kinematic range covered by the data input to the QCD 
fit is 3 x 10~ 5 < x < 0.7 and 1< Q 2 < 5000 GeV 2 . 

The QCD predictions for the F% structure functions are obtained by solving the DGLAP 
evolution equations |j] at NLO in the MS scheme [52]. These equations yield the quark 
and gluon momentum distributions (and thus the structure functions) at all values of 
Q 2 provided they are given as functions of x at some input scale Qq. In this analysis we 
adopt the so-called fixed flavour number scheme where only three light flavours (u, d, s) 
contribute to the quark density in the proton. In this scheme the assumption is made 
that the charm and bottom quarks are produced in the hard scattering process and 
the corresponding structure functions F% and F\ are calculated from the photon-gluon 
fusion process including NLO corrections |5l 



As will be explained later, the input scale is chosen to be Q\ = 7 GeV 2 . The gluon 
distribution (xg), the sea quark distribution (xS) and the difference of up and down 
quarks in the proton (xA u d) are parameterised as 

xg(x,Q 2 ) = A g x^{l-x)^{l + lg x) 

xS(x,Ql) = 2x(u + d+s) = A s X Ss (l-xf s (l + £ s y/x~ + ls x) (7) 
xA ud (x,Ql) = x{u + u) - x(d + d) = A a x 5a (1 - x) VA . 

The input valence distributions xu v = x(u — u) and xd v = x(d — d) at Qq are taken 
from the parton distribution set MRS(R2) |pi| . As for MRS(R2) we assume that 
the strange quark distribution is a given fraction K s = 0.2 of the sea at the scale 

6 Data from E665 are not included in the fit. They are important at low x and Q 2 < f GeV 2 but 
of much lower statistical weight at larger x compared to BCDMS and NMC. We have checked that 
including E665 data within the cuts described does not change the nominal fit result. 
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Q 2 = 1 GeV 2 . The gluon normalisation, A g , is fixed by the momentum sum rule, 

/ (xg + xS + xu v + xd v )dx = 1. There are thus 11 free parameters in the fit. 
Jo 

The input value for the strong coupling constant is set to a s (M§) = 0.118 [55|. With a 
charm (bottom) threshold of Q c (b) = 1-5 (5) GeV this corresponds to values of the QCD 
scale parameter A^jg = (404, 343, 243) MeV for / = (3, 4, 5) flavours. In the calculation 
of the charm structure function F% the charm mass is taken to be m c = 1.5 GeV; the 
contribution from bottom is found to be negligible in the kinematic range covered by 
the data. In the QCD evolutions and the evaluation of the structure functions the 
renormalisation scale and the mass factorisation scale are both set equal to Q 2 . 
The QCD evolutions and the structure function calculations are done with the program 
QCDNUM |56[ |. The QCD evolution equations are written in terms of quark flavour 
singlet and non-singlet and gluon momentum distributions. The quark non-singlet evo- 
lution is independent of the gluon. The quark singlet distribution is defined as the sum 
over all quark and anti-quark distributions 

= [xqi(x) + xqi(x)] (8) 

i=u,d,s 

and its evolution in Q 2 is coupled to that of the gluon distribution. At small values of x, 
xT, is dominated by the contribution from the qq sea xS. Note that for data with Q 2 < 
Qq = 7 GeV 2 , backwards evolutions in Q 2 are performed. The x 2 minimisation and the 



calculation of the covariance matrices are based on MINUIT [57]. In the definition of 
the x 2 only statistical errors are included and the relative normalisations of the data 
sets are fixed at unity. 

The fit yields a good description of the data as shown in Fig. 10 where we plot the Q 2 
dependence of the proton structure function in the x range covered by the ZEUS94 
data. The characteristic pattern of scaling violations can be seen clearly from this plot, 
with F2 at low values of x rising as Q 2 increases. The quality of the fit to ZEUS data 
at low Q 2 is also shown by the full line in Fig. ||. Adding the statistical and systematic 
errors in quadrature gives a x 2 of 1474 for 1120 data points and 11 free parameters. 
We have also checked that the gluon obtained from this fit to scaling violations gives 
values of F% in agreement with the ZEUS measurements [58]. The values of the fitted 
parameters are given in Table || of the Appendix. 

Fig. [ll| shows the gluon momentum distribution as a function of x for Q 2 at 1, 7 and 

20 GeV 2 . The following sources contribute to the inner shaded error bands displayed in 

the figure (for each source we give in brackets the relative error Ag/g at x = 5 x 10 -5 , 
Q2 = 7 GeV 2) : 



1. The statistical error on the data 

2. The experimental systematic errors, (13%), which are propagated to the covariance 



matrix of the fitted parameters using the technique described in [59|. In total 26 
independent sources of systematic error are included. Proper account is taken 
of the correlations between the systematic errors of the NMC datasets and for 



BCDMS the procedure of [ 50 1 is followed. Normalisation errors of all data sets are 
also included. 
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3. The uncertainties on the strong coupling constant Aa s (M§) = 0.005 (8%), on 
the strange quark content of the proton AK S = 0.1 (1%) and on the charm mass 
Am c = 0.2 GeV (1%). 

Adding errors from (1), (2) and (3) together in quadrature gives a total contribution of 
16%. [] 

In addition to the above sources of error a 'parameterisation error' (10%) is obtained by 
repeating the fit with: 

4. The addition of statistical and systematic errors in quadrature in the definition of 
the x 2 instead of taking statistical errors only. 

5. The input scale set to Q% = 1 and 4 GeV 2 instead of 7 GeV 2 . 

6. An alternative parameterisation of the gluon density: 

3 



xg(x,Q ) = A g (l 



l + C n T n (y) 



n=i 



0) 



where T n (y) is a Chebycheff polynomial of the first kind |30[ and y = alnx + b 
with the coefficients a, b adjusted such that x £ [10 -6 , 1] maps onto y S [—1,1]. 
This parameterisation is flexible enough to describe the rapid change with Q 2 
(see Fig. [TT|) of the shape of the gluon density. Furthermore, Eq. (||) allows the 
gluon density to become negative at low x whereas Eq. (0) imposes the constraint 
xg{x) > as x — > 0. 

Taking all combinations of the alternatives described in 4., 5. and 6. above in addition to 
the nominal settings, twelve fits are performed and the parameterisation error is defined 
as the envelope of the resulting set of quark and gluon distributions. All these fits yield 
similar values of x 2 - The nominal fit (ZEUSQCD) is taken to be that which gives a 
curve that is roughly at the centre of the error bands for all x and Q 2 . This also defines 
the choice of Qq = 7 GeV 2 . The outer hatched error bands in Fig. 11 correspond to the 
total error now including the parameterisation error added in quadrature with the other 
errors. At x = 5 x 10~ 5 , Q 2 = 7 GeV 2 the total Ag/g = 19%. 



The three left-hand plots of Fig. |12| show the distributions for xE and xg as functions 
of x for Q 2 at 1, 7 and 20 GeV 2 . The error bands shown correspond to the quadratic 
sum of all error sources. It can be seen that even at the smallest Q 2 the quark singlet 
distribution, xT,, is rising at small x whereas the gluon distribution, xg, has become 
almost flat, indeed compatible with zero. This behaviour has also been found by others, 
for instance Martin et al. (MRST) |4(| in their recent global determination of parton 
densities. At Q 2 = 1 GeV 2 the gluon distribution is poorly determined and can, within 
errors, be negative at low x. In the simplest form of the parton model (and leading 
order QCD) this would clearly be unphysical and while it is known that at NLO in the 
MS scheme a positive parton density will remain positive for forward evolution in Q 2 
there is no such constraint for backwards evolution . A negative gluon distribution is 
therefore not necessarily in contradiction with perturbative NLO QCD provided cross- 
sections or structure functions calculated from the parton distributions are positive for 

7 This combination of errors ('HERA standard errors') is often used by the HI and ZEUS experiments 
when discussing xg. 
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all x and Q 2 in the fitted kinematic domain. We have verified that this is the case for 
F$ and the longitudinal structure function Fl$ 

7 The transition region and NLO QCD 
7.1 The NLO QCD fit at low Q 2 

It is now widely observed that NLO DGLAP QCD fits give good descriptions of F2 data 
down to Q 2 values in the range 1 — 2 GeV 2 . For such fits to be valid one assumes: 

1. the validity of the DGLAP QCD formalism; 

2. NLO is sufficient even though as is becoming large (e.g. for this analysis as = 0.46 
at 1 GeV 2 ); 

3. that no higher twist terms, shadowing or other non-perturbative effects contribute 
to F 2 . 

In this paper we have also deliberately made the minimum of assumptions about the 
low x functional form of the parton distributions at Qq. We require only that they must 
tend to zero as x — ► 1 and that the flavour and momentum sum-rules are respected. 
To investigate the stability of our results at low Q 2 we have repeated the full QCD fit 
and error evaluation procedure on the same data as in Sec. ^ but with the minimum Q 2 
cut {Q 2 m in) raised to 4 GeV 2 . The quality of the fit is much as before, x 2 of 1242 for 943 
data points and 11 parameters (statistical and systematic errors added in quadrature). 
The resulting xS and xg parton distributions are shown in the three right-hand plots 
of Fig. [l^. Qualitatively the features shown by the standard fit (left-hand plots) are 
unchanged. The rising xT, distribution at low x remains and the sea dominates the 
gluon at small x and the lowest Q 2 value. In more detail: 

• The parton densities from the central fits with Q 2 ^ = 1 and 4 GeV 2 are very 
similar. 

• Except at the lowest Q 2 , the precision of the determination of is not much 
reduced. Even at Q 2 = 1 GeV 2 , xT, is reasonably well determined for x > 3 x 10 -4 , 
at smaller x values there are insufficient data to constrain the fit. 

• At all Q 2 values shown in the right-hand plots the precision of the determination 
of xg for x < 10 -3 is worse. 

The increase in the error bands for the gluon density when the Q 2 ^ cut is increased 
shows the importance of the SVX95 data at low x and low Q 2 in determining xg in this 
region. 

To investigate if there is a technical lower limit to the NLO QCD fit (in the sense that 
the fit fails to converge or gives a very bad x 2 )> we extend the QCD fit into the region 
covered by the ZEUS BPC95 data by lowering the Q 2 min cut to 0.4 GeV 2 . The fit gives 
an acceptable description of the data with F% positive and Fl (calculated as in fl62| ) only 
slightly negative at x ~ 5 x 10~ 3 and Q 2 = 0.5 GeV 2 . We therefore conclude that we do 

8 Here and in the following Fl is calculated using the QCD calculation of Altarelli and Martinelli |i2j . 
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not observe a significant breakdown of the NLO DGLAP description in the kinematic 
range explored. However, we stress that, within the present experimental accuracy, F2 
data by itself cannot validate the assumptions noted at the beginning of this section 
and that other information such as precise measurements of F% in DIS, measurements 
of other hard processes or more theoretical input is required. f\ 

7.2 F 2 slopes and models 

In order to clarify what we have learned about the transition region from the ZEUS 
Regge fit and the NLO DGLAP fits we return to the F2 slopes. 

Fig. 0, showing the ZEUS+E665 A e // data of Sec. 5J, also shows the calculation from 



the ZEUSREGGE fit of Sec. | and for Q 2 > 1 GeV 2 the result from the ZEUSQCD 
fit of Sec. ^ In the BPC region, the ZEUSREGGE (full line) calculation gives A e // 
somewhat higher than that given by the original DL Regge fit (dashed line), but it 
is still below the data. This is largely because the ZEUSREGGE fit includes low W 
photoproduction data and thus gives a lower ctp than the value of 1.141 from the Regge 
fit to the extrapolated ZEUS BPC ctq P data alone. At larger Q 2 values, ZEUSQCD 
(full line) gives a good account of the trend and normalisation of X e jf while GRV94 
(dashed line) tends to predict a larger value of A e jj, which means a steeper rise of F2 
as x decreases, than that determined from the data. 

Fig. I shows the ZEUS dF 2 /dhiQ 2 data together with the same two calculations, 
ZEUSREGGE at very low x and Q 2 values and ZEUSQCD for Q 2 > 1 GeV 2 . The 
ZEUSREGGE points show much the same trend as that of the original DL model. At 
Q 2 values between 1 and 5 GeV 2 , the ZEUSQCD points are qualitatively different from 
those of GRV94. The ZEUSQCD values now follow the 'turn over' of the slope data 
around x ~ 3 x 10~ 4 . This has also been found by MRST in 46] where they compare 



their global fit to (IF2/ 'd\nQ 2 data from HI, ZEUS and the NMC experiments. The rise 
in F2 at low x and the 'turn over' in dF2 /d In Q 2 reflect the different behaviours of the 
gluon and qq sea distributions at Q 2 ~ 1 GeV 2 . For very small x the sea continues to 
rise (5 S < in the notation of Eq. (|7|)) whereas the gluon rises significantly less steeply 
or even tends to zero (5 g > 0). 

In contrast, for the GRV94 parton densities GRV assume that at their low starting 
scale of 0.34 GeV 2 both the gluon and the sea distributions are non-singular. All the 
rise of F2 at low x for Q 2 > 1 GeV 2 is then generated through the DGLAP evolution 
equations. In a recent paper GRV |64j have revisited their 'dynamical parton model' in 
the light of the HERA 1994 data and have produced a new set of parton distributions 
- GRV98. They find that they can correct most of the discrepancies between GRV94 
and the HERA data by a slight increase in the starting scale from 0.34 to 0.4 GeV 2 and 
by using a lower value of as(M|) = 0.114 rather than the value of 0.118 used in the 
ZEUS NLO QCD fit. They acknowledge that if they use a larger value of as then their 
starting scale has to be increased to around 1 GeV 2 and they have to accept a rising qq 
sea distribution at the starting scale. 



9 Some of the features that we find, such as the suppression of the gluon density at small Q 2 , and the 
ability of the NLO DGLAP formalism to fit very low Q 2 data have been noted previously by Lopez et 
al. @. 
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8 Summary and conclusions 



In this paper we have presented the measurement by ZEUS of F2 in the Q 2 region 
0.6 - 17 GeV 2 (SVX95), which fills the gap between the very low Q 2 BPC95 data 
(0.11 - 0.65 GeV 2 ) and the large 1994 data sample (3.5 - 5000 GeV 2 ). We have shown 
that the BPC data may be described successfully by non-perturbative approaches: a 
simple generalised vector dominance model for the Q 2 dependence and a two component 
Regge model for the W 2 dependence. For Q 2 > 0.9 GeV 2 these approaches fail to 
describe the dominant feature of the data, which is the rapid rise of F2 at small x. 
We have studied the transition region by fitting F2 = Ax" e ff using ZEUS and E665 
data with x < 0.01 in the Q 2 range 0.15 - 250 GeV 2 . For Q 2 > 0.9 GeV 2 the data are 
not compatible with the Q 2 independence of A e yj, as expected for a dominant Pomeron 
term in conventional Regge theory, but are well described by the ZEUS NLO QCD fit. 
The slope (IF2/ d\xiQ 2 has been calculated from ZEUS F% data in the range 2 x 10~ 6 < 
x < 0.2. Assuming pQCD, F2 at low x is largely determined by the qq sea density 
whereas dF2/d\nQ 2 is given by the gluon density. As x decreases the slope values 
increase until at x ~ 3 x 10 -4 there is a turn over and for smaller x values the slope 
values decrease. 

To study the behaviour of the parton densities in more detail we have performed a NLO 
DGLAP QCD fit to ZEUS and fixed target data with Q 2 > 1 GeV 2 and 3 x 10~ 5 < x < 
0.7. A good description of the F2 data over the whole range of Q 2 from 1 to 5000 GeV 2 
is obtained. Around the lower Q 2 limit of the fit we find that the qq sea distribution 
is still rising at small x, whereas the gluon distribution is strongly suppressed. These 
findings are incompatible with the hypothesis that the rapid rise in F% is driven by the 
rapid increase in gluon density at small x from parton splitting alone. These features 
remain true if the lowest Q 2 for which data is included in the fit is raised from 1 to 
4 GeV 2 . 

From the ZEUS QCD fit we also obtain a much improved determination of the gluon 
momentum density compared to the previous determination by ZEUS ||. Full account 
has been taken of correlated experimental systematic errors and the uncertainty in the 
form of the input gluon distribution function has also been estimated. At Q 2 = 20 GeV 2 
and x = 5 x 10 -5 the total fractional error on the gluon density has been reduced from 
40% to 10%. 

We have used NLO pQCD and the simplest non-perturbative models to study the tran- 
sition region in Q 2 . We find, for Q 2 > 0.9 GeV 2 , that the wholly non-perturbative 
description fails. Although pQCD may not be valid at such low scales, we have not been 
able to find a lower limit at which the NLO DGLAP fit breaks down conclusively. 
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Appendix 



The following tables are also available via the ZEUS collaboration home page, 
http:/ / www-zeus.desy.de/ , Fortran routines and data files to calculate the parton dis- 
tributions from the ZEUS NLO QCD fits are also available from this site. 
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bin 


X 


Q 2 ( GeV 2 ) 


F 2 


error stat. 


1 


1 


2 


• 10 


-5 


0.6 





531 


±0.0301^ 


2 


1 


9 


• 10 


-5 


0.9 





653 


±0.030l"g 


3 


2 


6 


• 10 


-5 


0.9 





655 


±0.02318 


4 


3 


6 


• 10 


-5 


0.9 





574 


±0.02218 


5 


2 


8 


• 10 


-5 


1.3 





716 


±0.02918 


6 


3 


8 


• 10 


-5 


1.3 





717 


±0.02318 


7 


5 


3 


• 10 


-5 


1.3 





715 


±0.02218 


8 


8 





• 10" 


-5 


1.3 





626 


±0.01718 


9 


1 


6 


• 10 


-4 


1.3 





559 


±0.01618 


10 


3 


9 


• 10 


-5 


1.9 





915 


±0.03618 


11 


5 


3 


• 10 


-5 


1.9 





852 


±0.028l[] 


12 


7 


5 


• 10 


-5 


1.9 





742 


±0.025l|] 


13 


1 


1 


• 10 


-4 


1.9 





740 


±0.02318 


14 


2 


2 


• 10 


-4 


1.9 





636 


±0.01918 


15 


5 


4 


• 10 


-5 


2.5 





964 


±0.03918 


16 


7 


3 


• 10 


-5 


2.5 





909 


±0.02918 


17 


1 





• 10 


-4 


2.5 





870 


±0.02718 


18 


1 


6 


• 10 


-4 


2.5 





771 


±0.02318 


19 


3 





• 10" 


-4 


2.5 





725 


±0.02318 


20 


7 


3 


• 10 


-5 


3.5 


1 


140 


±0.05018 


21 


9 


8 


• 10 


-5 


3.5 





995 


±0.03518 


22 


1 


4 


• 10 


-4 


3.5 





945 


±0.03118 


23 


2 


1 


• 10 


-4 


3.5 





839 


±0.02618 


24 


4 


1 


• 10 


-4 


3.5 





663 


±0.02118 


25 


1 





• 10 


-4 


4.5 


1 


160 


±0.04318 


26 


1 


7 


• 10 


-4 


4.5 


1 


005 


±0.03118 


27 


2 


8 


• 10 


-4 


4.5 





912 


±0.03018 


28 


5 


4 


• 10 


-4 


4.5 





664 


±0.02218 


29 


1 


4 


• 10 


-4 


6 


1 


305 


±0.05518 


30 


2 


2 


• 10" 


-4 


6 


1 


086 


±0.03918 


31 


3 


7 


• 10 


-4 


6 


1 


002 


±0.03818 


32 


7 


1 


• 10 


-4 


6 





753 


±0.02818 


33 


1 


8 


• 10 


-4 


7.5 


1 


210 


±0.05818 


34 


3 


6 


• 10 


-4 


7.5 





918 


±0.03618 


35 


8 


8 


• 10 


-4 


7.5 





810 


±0.04018 


36 


2 


2 


• 10 


-4 


9 


1 


296 


±0.06ll8 


37 


4 


3 


• 10 


-4 


9 





957 


±0.03818 


38 


1 


1 


• 10 


-3 


9 





818 


±0.03918 


39 


2 


8 


• 10 


-4 


12 


1 


306 


±0.05718 


40 


5 


6 


• 10 


-4 


12 


1 


084 


±0.04018 


41 


1 


4 


• 10" 


-3 


12 





956 


±0.04418 


42 


3 


9 


• 10 


-4 


17 


1 


311 


±0.05718 


43 


7 


9 


• 10 


-4 


17 


1 


087 


±0.03818 


44 


1 


9 


• 10 


-3 


17 





931 


±0.04518 



-sys 



-0.070 



-0.044 



-0.023 
H0.034 
-0.022 
h0.019 
-0.028 



-0.079 



-0.074 



-0.038 
1-0.043 
-0.061 



h0.051 
-0.041 
h0.046 
-0.077 



F 2 (F L = 0) 
0.517 
0.632 
0.644 
0.570 
0.687 
0.703 
0.708 
0.624 
0.558 
0.873 
0.833 
0.734 
0.736 
0.635 
0.916 
0.887 
0.860 
0.767 
0.724 
1.081 
0.970 
0.933 
0.834 
0.661 
1.108 
0.990 
0.907 
0.663 
1.247 
1.070 
0.996 
0.751 
1.167 
0.910 
0.808 
1.251 
0.948 
0.816 
1.265 
1.073 
0.954 
1.270 
1.077 
0.927 



Table 1: Values of x, Q 2 , F 2 , statistical and 
shifted vertex analysis (SVX95). 18 



systematic errors from the ZEUS 1995 



bin 


W 


e 


< 


stat. err. 


sys. err. 




[GeV] 






[H 


[H 


1 


104 


0.99 


156.2 


±5.3 


±16.1 


2 


134 


0.98 


166.1 


±5.2 


±11.0 


3 


153 


0.96 


174.7 


±4.9 


±12.9 


4 


173 


0.92 


175.5 


±5.0 


±11.7 


5 


190 


0.88 


181.8 


±4.7 


±12.8 


6 


212 


0.80 


186.8 


±4.8 


±13.5 


7 


233 


0.69 


192.5 


±4.7 


±13.3 


8 


251 


0.55 


204.8 


±5.6 


±17.0 



Table 2: Values of a^ p together with statistical and systematic errors from the GVMD 
extrapolation, with ox = 0, of ZEUS BPC95 data. W is the 7*p centre-of-mass energy 
and e is the polarisation of the virtual photon. 
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1 ■ 

bin 


Q z ( GeV z ) 


X 


rain 




X 


max 




( 


x) 




X eff 


error stat. 


1 


0.15 


2.4 


10" 


-6 


4.2 


io- 


-71 

-b 


3.2 


io- 


—p 

-D 


0.134 


1 TTT 

±0.0771^ 


2 


0.2 


3.2 


10" 


-6 


8.5 


io- 


-6 


5.5 


io- 


-6 


0.128 


±0.036±^ 


3 


0.25 


4.6 


io- 


-6 


1.8 


io- 


-3 


1.3 


10" 


-5 


0.136 


±0.016lg 


4 


0.3 


6.7 


io- 


-6 


2.5 


io- 


-3 


1.1 


10" 


-4 


0.110 


±0.007lg 


5 


0.4 


1.1 


10" 


-5 


5.2 


io- 


-3 


6.8 


io- 


-4 


0.115 


±0.005i£ 


6 


0.5 


2.1 


io- 


-5 


4.6 


10" 


-5 


2.9 


10" 


-5 


0.193 


±0.050±g 


7 


0.65 


1.2 


io- 


-5 


5.2 


io- 


-3 


1.8 


io- 


-3 


0.114 


±0.007ig 


8 


0.9 


1.9 


io- 


-5 


8.9 


io- 


-3 


3.8 


10" 


-3 


0.147 


i n 

±0.005lg 


9 


1.3 


2.8 


10" 


-5 


8.9 


io- 


-3 


2.4 


io- 


-3 


0.152 


±0.006l[J 


10 


1.9 


3.9 


10" 


-5 


8.9 


io- 


-3 


1.8 


io- 


-3 


0.160 


i n 

±0.008i[J 


11 


2.5 


5.4 


io- 


-5 


8.9 


io- 


-3 


1.3 


10" 


-3 


0.179 


±0.010^ 


12 


3.5 


6.3 


10" 


-5 


8.9 


io- 


-3 


6.1 


io- 


-4 


0.178 


±0.012ig 


13 


4.5 


1.0 


io- 


-4 


5.4 


10" 


-4 


2.0 


10" 


-4 


0.261 


±0.020l[J 


14 


6.5 


1.0 


io- 


-4 


4.0 


io- 


-3 


7.9 


10" 


-4 


0.192 


±0.005ig 


15 


8.5 


1.6 


10" 


-4 


1.6 


io- 


-3 


6.7 


10" 


-4 


0.261 


±0.0121° 


16 


10 


1.6 


10" 


-4 


6.3 


io- 


-3 


1.6 


io- 


-3 


0.226 


±0.007lg 


17 


12 


2.5 


10" 


-4 


2.5 


io- 


-3 


9.7 


io- 


-4 


0.250 


i n 

±0.014jlg 


18 


15 


2.5 


io- 


-4 


6.3 


10" 


-3 


1.9 


10" 


-3 


0.249 


iO.OlOlg 


19 


18 


3.9 


io- 


-4 


6.3 


io- 


-3 


2.1 


io- 


-3 


0.274 


±0.010lg 


20 


22 


4.0 


10" 


-4 


1.0 


io- 


_2 


3.4 


io- 


Q 

o 


0.273 


±0.0101q 


21 


27 


6.3 


10" 


-4 


6.3 


io- 


-3 


2.7 


io- 


-3 


0.332 


±0.015l[| 


22 


35 


6.3 


io- 


-4 


1.0 


io- 


-2 


3.7 


io- 


-3 


0.308 


±0.01618 


23 


45 


1.0 


10" 


-3 


1.0 


io- 


-2 


3.7 


10" 


-3 


0.289 


±0.02018 


24 


60 


1.0 


io- 


-3 


1.0 


io- 


-2 


4.7 


io- 


-3 


0.342 


±0.02118 


25 


70 


1.6 


io- 


-3 


1.0 


10" 


-2 


4.8 


10" 


-3 


0.306 


±0.03218 


26 


90 


1.6 


10" 


-3 


1.0 


io- 


-2 


5.3 


io- 


-3 


0.273 


±0.04018 


27 


120 


2.5 


10" 


-3 


1.0 


io- 


-2 


6.2 


10" 


-3 


0.380 


±0.05918 


28 


150 


2.5 


10" 


-3 


1.0 


io- 


-2 


6.3 


io- 


-3 


0.258 


±0.07918 


29 


200 


4.0 


io- 


-3 


1.0 


10" 


-2 


7.0 


10" 


-3 


0.250 


±0.12618 


30 


250 


4.0 


10" 


-3 


1.0 


10" 


-2 


7.9 


10" 


-3 


0.452 


±0.18418 



+sys 
-sys 



-0.083 
-0.040 



-0.019 



-0.018 
-0.042 
-0.041 



-0.020 
-0.016 
-0.014 



-0.019 



-0.009 
-0.017 
-0.018 
-0.023 
-0.046 
-0.027 



-0.003 
-0.008 
-0.000 



Table 3: Values of the slope X e ff = din F2 / d\n(l / x) and their errors, calculated from 
fitting F 2 = Ax~ x <>ff at fixed Q 2 to ZEUS and E665 data with x < 0.01. The columns 
labelled x m i n , x max and (x) give the minimum, maximum and average values of a; in a 
bin. (x) is calculated as described in Sec. [O] of the text. 
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bin 
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Qmin \ tjeV ) 




/r> 2 \ 
\Q ) 


at '2/ aiiUcj 


error stat. 


1 


2.1-10" 6 


0.11 


0.15 


0.12 


0.135 


±0.029lo 


2 


3.M0~ 6 


0.15 


0.2 


0.16 


0.198 


±0.028ig 


3 


4.6-10~ 6 


0.15 


0.25 


0.2 


0.174 


±0.017lo 


4 


7.3-10~ 6 


0.2 


0.3 


0.23 


0.191 


±0.0191q 


5 


1.2-10~ 5 


0.25 


0.6 


0.29 


0.265 


±0.017±q 


6 


2.0-10~ 5 


0.3 


0.9 


0.37 


0.297 


±0.0181q 


7 


3.3-10~ 5 


0.3 


1.9 


0.52 


0.312 


±0.009to 


8 


6.3-10 -5 


0.5 


3.5 


1.1 


0.365 


io.oioio 


9 


1.0-10" 4 


1.3 


6.5 


2.5 


0.379 


±0.0181q 


10 


1.6-10 -4 


1.3 


10 


3.8 


0.387 


±0.0131q 


11 


2.5-10~ 4 


1.9 


15 


5.2 


0.368 


±0.013lo 


12 


4.0-10" 4 


3.5 


22 


8.8 


0.429 


±0.0181q 


13 


6.3-10" 4 


4.5 


35 


10 


0.404 


±0.0151q 


14 


1.0-10~ 3 


6.5 


60 


13 


0.315 


±0.014lg 


15 


1.6-10~ 3 


6.5 


90 


14 


0.262 


±0.0121q 


16 


2.5-10~ 3 


6.5 


150 


19 


0.227 


±0.0091° 


17 


4.0-10~ 3 
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26 


0.139 


±0.008io 


18 


6.3-10 


10 


A cr n 

450 


O A 

24 
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±0.008_ 
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1.0-10~ 2 


22 


800 
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±0.012±g 


20 


1.6-10" 2 
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0.059 


±0.005ig 


21 


2.5-10" 2 


22 
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45 


0.061 


±0.008ig 


22 


4.0-10~ 2 
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2000 


167 


0.037 


±0.007ig 


23 


8.1-10~ 2 


10 


5000 


156 


0.018 


±0.004±8 


24 


0.2 


90 


5000 


388 


0.008 


iO.OlOtg 



+sys 
-sys 
-0.U30 



-0.022 
-0.032 
-0.031 



-0.040 



-0.020 
-0.015 



-0.011 



-0.010 
i-0.019 
-0.045 



Table 4: Values of the slope (IF2/ d\nQ 2 and their errors calculated from fitting ZEUS 
F2 data to the form a + 61nQ 2 . The columns labelled QmimQmax an d {Q 2 } give the 
minimum, maximum and average values of Q 2 in a bin. (Q 2 ) is calculated as described 



in Sec. 5.2 of the text. 



Parameter 


xg(x) 


xS(x) 


xA ud (x) 


A 


1.77 


0.520 


6.07 


5 


-0.225 


-0.241 


1.27 


77 


9.07 


8.60 


3.68 


e 




0.290 




7 


3.00 


8.27 





Table 5: Values for the parameters of the nominal ZEUS NLO QCD fit (ZEUSQCD) at 
the starting scale Qq = 7 GeV 2 . The parameters are defined in Eq. (^) of the text. 
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Figure 1: The (x, Q 2 ) plane showing the regions covered by the ZEUS data sets BPC95, 
ZEUS94 and SVX95 together with regions covered by the fixed target experiments E665, 
BCDMS and NMC. 
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Figure 2: Distributions from the SVX95 analysis showing the comparison between data 
and simulation: (a) 8 or E — Pz as defined in Sec. [D^; (b) the Z-position of the primary 
vertex; (c) the positron scattering angle 6 e ; (d) the energy, E' e , of the scattered positron; 
(e) log 1Q Q 2 ; (f) log 10 y- In all cases the data are represented by filled circles, the 
simulation by the open histograms and the photoproduction background (calculated 
from the MC described in Sec. 3.1) by the filled triangles. Both MC calculations are 
normalised to the luminosity of the data. 
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Figure 3: The ZEUS SVX95 F2 data as a function of x for different Q 2 bins together 
with previous ZEUS data (ZEUS 1994 §, ZEUS BPC 1995 §), data from HI SVX95 



rrj] and fixed target data (E665 [32]). Error bars correspond to the statistical and 
systematic errors added in quadrature. The overall normalisation errors are not shown. 
The curves shown are (dotted) the Donnachie-Landshoff Regge model |22j and (full) the 
ZEUS NLO QCD fit. 
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Figure 4: The total 7*p cross-section, <r^ ot p , as a function of W 2 at different Q 2 [GeV 2 ]. 
The data of this analysis (ZEUS SVX95) are shown together with previous ZEUS and 
HI data (ZEUS94, ZEUS BPC95, HI SVX95). Also the total photoproduction cross- 
sections from ZEUS and HI and from fixed target experiments at lower W are shown. 
Predictions from the DL [p2j | and GRV94 j|] models are indicated by the dotted and 
dashed curves, respectively. 
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Figure 5: Upper plot. ZEUS BPC measurements of the total cross-section o~t + eai 
in bins of W as a function of Q 2 and the GVMD fit to the data. Lower plot, crj^j. as 
a function of W 2 . The ZEUS BPC95 points are those from the GVMD extrapolation 
(<7q P ). Also shown are direct measurements of the total photoproduction cross-section 
from HI, ZEUS and earlier experiments at low energies. The curves show Regge fits: 
the original DL fit [^] to the low W data (dotted); the Pomeron only fit to the BPC 
<Tq P data (dashed) and the Pomeron+Reggeon fit to the low W and BPC <7q P data (full). 
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Figure 6: Low Q 2 F2 data for different Q 2 bins together with the ZEUS Regge fit 
(dashed curves) to the BPC95 data as described in Sec. Also shown at larger values 
of Q 2 is the ZEUS NLO QCD fit (full curves) as described in Sec. §. 
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Figure 7: A e // = d\n F2 / d\n(l / x) as a function of Q 2 calculated by fitting F2 = Ax~ e ff 
to ZEUS and E665 data with x < 0.01. The inner error bar shows the statistical error 
and the outer the total statistical and systematic error added in quadrature, (x) is 



calculated as described in Sec. 5JJ . The DL and GRV94 calculations, shown as points 
linked by dashed lines, are from the Donnachie-Landshoff Regge fit |22| and the GRV94 
NLO QCD fit §§, respectively. The ZEUSREGGE and ZEUSQCD calculations, shown 
as points linked by full lines, are from the ZEUS Regge and NLO QCD fits described in 
Sees |] and || of the text respectively. In all cases the points are obtained using the same 
weighted range of for the experimental data. 



31 



ZEUS 1995 



IN 
fa 



_ x=2.1E-06 

. <Q*>=0.12GeV 

^ • ZEUS BPC 1995 
■ ZEUS SVTX 1995 
A ZEUS 1 994 




_ x=3.1E-06 

. <0 2 >=0.16GeV J 

V 




H x=4.6E-06 
_ <Q 2 >=0.2GeV 2 




_ x=7.3E-06 

. <0 2 >=0.23GeV 2 




_ x=1.2E-05 

. <Q 2 >=0.29GeV 2 

V 

I , I , I , I 


- x=2E-05 

. <Q 2 >=0.37GeV 

■ 

V 

I I I , I , 


_ x=3.3E-05 

. <Q 2 >=0.52GeV 2 

/ 

, , , I , I 


_ x=6.3E-05 
. <Q 2 >=1.1GeV 2 

: k 

- / 

i , , i , i 


_ x-0.0001 

_ <Q 2 >=2.45GeV 2 

: / 

I 

, I , I , I , I 


_ x-0.0001 6 

. <tf> =3.8506^ 

/ 

, I , I , I , I 


. x-0.00025 

. <Q 2 >=5.24G^V 2 

M 

, I , I , I , I 


. x-0.0004 

. <Q 2 >=8.76GeV 2 

: / 

m 

, i , i , i , i 


. x=0.00063 

_ <Q*> = 10.63G<*" 

* 

, I , I , I , I 


. x=0.001 

_ <Q*> = 13.52Ge\£ 

: J 

L^j — u — ^_ 




_ x =0.00 16 

_ <Q 2 >=14.57GeV 2 

J> 

■ / 

, I , I , I , I 


_ x=0.0025 

_ <0 2 >=19.26GeV 2 

1 

■ / 

, i , i , i , i 


. x=0.004 

. <Q 2 >=26.15GeV 2 

I i I , I I 


_ x=0.0063 

. <Q 2 >=24.58GeV 2 ' 

- J* 

I I 


i 

I 


_ x=0.01 

. <Q 2 >=56.27GeV 2 

i , i i i i 


_ x=0.016 

_ <Q J >=19.37GeV 2 


. x=0.025 

_ <0 2 >=45.01GeV 2 


_ x=0.04 

_ <Q J > = 167.61GeV 2 




. x=0.081 

_ <Q 2 >=156.59GeV 2 


. x=0.2 

_ <0 2 >=388.03GeV 2 


, I , I , I , I 


, I , I 










i io 3 i io 3 i io 3 i io 3 

Q 2 (GeV 2 ) 



F2 as a function of Q 2 in bins of x from the ZEUS data sets BPC95, SVX95, 
The linear fits F2 = a + 61nQ 2 are also shown. The values of (Q 2 ) given in 



Figure 8: 
ZEUS94. 

the plots are calculated as described in Sec. 1^2 
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Figure 9: dF2/d In Q 2 as a function of x calculated by fitting ZEUS F2 data in bins 
of x to the form a + b In Q 2 . The inner error bar shows the statistical error and the 
outer the total statistical and systematic error added in quadrature. (Q 2 ) is calculated 



as described in Sec. 5.2. The DL and GRV94 calculations, shown as points linked by 



dashed lines, are from the Donnachie-Landshoff Regge fit []22| and the GRV94 NLO 
QCD fit @, respectively. The ZEUSREGGE and ZEUSQCD calculations, shown as 
points linked by full lines, are from the ZEUS Regge and NLO QCD fits described in 
Sees |] and || of the text respectively. In all cases the points are obtained using the same 
weighted range of Q 2 as for the experimental data. 
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Figure 10: The proton structure function F2 versus Q 2 at fixed values of x. Data are 
from the ZEUS94 and SVX95 analyses and from the NMC and BCDMS fixed target 
experiments. The solid lines correspond to the QCD fit described in the text. For clarity 
an amount Cj = 13.6 — 0.6i is added to F2 where i = 1 (18) for the lowest (highest) x 
value. 
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Figure 11: The gluon momentum distribution xg(x) as a function of x at fixed values 
of Q 2 = 1,7 and 20 GeV 2 from the ZEUS QCD fit. The inner shaded bands show the 
'HERA standard' errors of Sec. [| The outer hatched bands indicate the quadratic sum 
of the 'HERA standard' and the 'parameterisation' errors. 
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Figure 12: The quark singlet momentum distribution, xT, (shaded), and the gluon 
momentum distribution, xg(x) (hatched), as functions of x at fixed values of Q 2 = 1, 
7 and 20 GeV 2 from the ZEUS QCD fit. The error bands correspond to the quadratic 
sum of all error sources considered for each parton density. The three left-hand plots 
show the results from the standard fit of Sec. ^ including F2 data with Q 2 > 1 GeV 2 ; 
the three right-hand plots show the corresponding results from the fit described in Sec. 
for which the data must satisfy Q 2 > 4 GeV 2 . 
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